On the Energy Increase in Space-Collapse Models 
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A typical feature of spontaneous collapse models which aim at localizing wavefunctions in space 
is the violation of the principle of energy conservation. In the models proposed in the literature 
the stochastic field which is responsible for the localization mechanism causes the momentum to 
behave like a Brownian motion, whose larger and larger fluctuations show up as a steady increase 
of the energy of the system. In spite of the fact that, in all situations, such an increase is small 
and practically undetectable, it is an undesirable feature that the energy of physical systems is 
not conserved but increases constantly in time, diverging for t oo. In this paper we show that 
this property of collapse models can be modified: we propose a model of spontaneous wavefunction 
collapse sharing all most important features of usual models but such that the energy of isolated 
systems reaches an asymptotic finite value instead of increasing with a steady rate. 
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I. INTRODUCTION 

As it is well-known, space-collapse models |1|, 
H H 11 H ES [H El [II aim at a solution of the 
macro-objectification or measurement problem in quan- 
tum mechanics by suitably modifying the Schrodinger 
equation with non-linear stochastic terms. One of the 
characteristic features of these models is the violation of 
energy-conservation for isolated systems; such a viola- 
tion is determined by the stochastic process responsible 
for the localization mechanism, which induces larger and 
larger fluctuations of the wavefunction in the momentum 
space 0: these increasing fluctuations, in turn, deter- 
mine the increase of the energy of the system [l5|. For 
typical values of the parameters, such an increase is very 
small and undetectable with present-day technology [l| ; 
still, one would wish to restore the principle of energy 
conservation within space-collapse models. 

In this paper we make one step towards this direction: 
we analyze a model of wavefunction space-collapse for 
which the energy of isolated systems does not increase 
indefinitely, but reaches an asymptotic finite value. An 
analogy with the quantum Brownian motion will show 
that the stochastic process acts like a dissipative medium 
which thcrmalizes the system to a fixed temperature (the 
temperature of the medium) and will suggest how to re- 



store perfect energy conservation. 
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The paper is organized as follows: after a brief re- 
view of the main features of dynamical reduction models 
(Sec. we introduce the collapse-model which is the 
subject of the paper fSec. Illljl . In Sec. IIVI we study the 
master equation for the statistical operator originating 
from the stochastic dynamics: this will provide the ra- 
tionale for the choice of the localization operator which 
defines the model. In Secs.lVlto lVIIII we will study in de- 
tail the most relevant properties of the model: we will 
analyze the time evolution of Gaussian wavefunctions 
(Sec. 01; the collapse mechanism and collapse probabil- 
ity fSec. IVI|) : we will see that the physical predictions of 
the model agree with very high accuracy with standard 
quantum mechanical predictions and, at the same time, 
the model reproduces classical mechanics at the macro- 
scopic level fSec. IVlfl . We will finally discuss the issue 
related to energy non conservation (Sec. IVIII(I and con- 
clude with some final remarks fSec. IIX|) . 
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II. STRUCTURE OF DYNAMICAL 
REDUCTION MODELS 

The typical structure of the evolution equation of col- 
lapse models is 1 : 



dipt 



Hdt + Vx (A-r t )dW t 



2r t A + r z t )dt 



with: 



(i) 



(2) 



The operator H is related to the standard quantum 
Hamiltonian of the system, while A is the reduction op- 
erator, i.e., the operator on whose cigcnmanifolds one 
wants to reduce the statevector, as a consequence of 
the collapse mechanism; the positive constant A sets the 
strength of the collapse mechanism. The stochastic dy- 
namics is governed by a standard Wiener process Wt de- 
fined on a probability space (il,J-,¥). Note that the 
equation is non-linear but preserves the norm of the stat- 
evector. 

In the literature on collapse models, the operator A is 
usually assumed to be self-adjoint; in such a case, and if 
one further assumes that it has only a discrete spectrum, 
it can be proven [3] that the form of the second and third 
term of Eq. which modify the standard Schrodinger 
evolution, is such that: 

1. The statevector collapses with respect to the "pre- 
ferred basis" generated by the operator A, i.e., for almost 
all realizations of the stochastic process there exists an 
eigenstate \a n ) of A (depending of course on the realiza- 
tion of the stochastic process) such that: 



for t 



(3) 



2. The average E[{P) t ] = E[(i/j t \P\ipt}] of any oper- 
ator P which commutes with A is constant in time, 
i.e., E[(^ t |P|-0 t )] = (ipo\P\tpo) . In particular, if P is a 
projection operator relative to an eigenmanifold of A, 
this together with property © implies that the proba- 
bility for the statevector to be reduced into the eigen- 
manifold associated to P is equal to (tpo\P\ipo), i.e., to 
the probability that standard quantum mechanics asso- 
ciates to the collapse, as a result of a measurement of 
the operator A. This is due to the fact that (ipt\P\ipt) 
turns out to be a martingale, thanks to the particular 
structure of Eq.Q, so that by the martingale property 
E[(^|P|^ t )] = (VoWo) 1- 



1 Of course, this form can be generalized, e.g. by adding a finite 
(or countable) number of operators Ai, each of which is coupled 
to a Wiener process Wi. Moreover, the Wiener processes may 
be complex instead of real, as assumed here. 



It is important to keep in mind that the above results 
are valid only when the standard Hamiltonian H either 
commutes with A or is equal to zero; in all other cases, 
such results are only approximate, the approximation de- 
pending on the value of A. 



III. THE MODEL 

In the literature, A has been mainly taken equal to the 
position operator q, or a function of q like in the continu- 
ous version Q of the original GRW model 1] , the reason 
being that the operator q is the most natural candidate 
for localizing wavefunctions in space. As anticipated in 
the previous section, one consequence of such a choice is 
that the energy of the system increases in time, diverg- 
ing for time going to infinity; it is then natural to wonder 
whether a different choice for A can preserve all most im- 
portant features of collapse models, but at the same time 
cure this energy non-conservation. This problem finds a 
partly positive solution by making the following Ansatz 2 
for A: 



A 



l H P > 



(4) 



where p is the momentum operator. Moreover, we define 
the operator H as follows: 



H 



Ho 



Xa 



{q,p}, 



(5) 



where Hq is the standard quantum Hamiltonian. In the 
following sections we will analyze the most relevant phys- 
ical properties of the model and we will focus our atten- 
tion to the case of a free particle: Hq — p 2 /2m, where m 
is the mass of the particle. 

The model is defined in terms of the two constants A 
and a; for reasons which will be clear in the following, we 
will assume them to vary with the mass of the particle 
as follows: 



A 



m 



A 



m 



(6) 



2 A localizati on o perator involving q and p has also been consid- 
ered in ref. Il (1 (note however that the form of the localization 
operator is different from ours) but with a different aim, i.e., that 
of studying whether the presence of a p-term instead of only a q- 
term can improve the localization mechanism. The authors prove 
that, for any physically reasonable choice of the parameters of 
their model, such term does not affect in an appreciable way the 
collapse mechanism. Here we show that a p-term is important as 
it can change the time evolution of the mean ener gy, a voiding it 
to increase constantly in time. The authors of ref. |l4l analyze a 
stochastic differential equation similar to our Eq. Q where both 
a q and a p term are present: they mainly focus their attention 
on the application of the formalism to the theory of open quan- 
tum systems and decoherent histories. One of their main result 
is a localization theorem which we will apply to our model to 
prove the collapse of wavefunctions to localized states. 
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where too is a reference mass which we choose to be equal 
to that of a nucleoli while Ao and ag are fixed constants 
which we take equal to: 



A 
on 



10 



-2 m -2 



- 18 m 2 . 



(7) 
(8) 



As it will be shown in Sec. IVIII this numerical choice 
for the parameters guarantees that the model reproduces 
almost exactly the physical predictions of standard quan- 
tum mechanics at the microscopic level and reproduces 
classical mechanics at the macroscopic level. 

Before concluding this section, we note that, in order to 
find the solutions of Eq. Q and to study their properties, 
it is convenient to consider also a linearized version of 
Eq. CD 40: 



d(j>t{x) 



i A 
-j-Hdt + VXAd^t - -A^Adt 



where £t is a standard Wiener processes defined on a new 
probability space (O, Q). In ref. 4] the relation be- 
tween the probability measures Q and F and the relation 
between the stochastic processes Wt and £t are discussed. 
Here we simply recall how one can use the above linear 
equation to find a solution of Eq. : 

of Eq. 10, with the initial 



1. Find the solution 
condition 4>q =ipo- 

2. Normalize the solution 3 : cj)t 

3. Make the substitution: dt; t 



+ ipt = <h/\\<k\\- 
dW t = d£ t - 2V\r t ; 



the wavefunction tp t thus obtained is a solution of Eq. 



IV. THE MASTER-EQUATION FOR THE 
STATISTICAL OPERATOR 

In order to better understand the modifications to the 
Schrodinger dynamics induced by Eq. Q and the mo- 
tivations for the precise choice of its form, apart from 
the martingale structure, and in particular in order to 
see why the choice (@J for A can partially cure the 
problem of the energy increase, it is worthwhile con- 
sidering the related equation for the statistical operator 
pt = E[\ipt)(tpt\], which is given by 



d 

dt Pt 



-J. [H ,pt] - - [q, [q,Pt]} 



Xa 2 
2fr? 



\P, \p,pt]\ - ll>vP>Pt}]> ( 10 ) 



that is the typical structure of master-equation for the 
quantum description of Brownian motion, where both 



friction and diffusion are taken into account and posi- 
tivity of the statistical operator is granted at all times. 
The obvious difference between Eq. I|10|) and the master- 
equation for quantum Brownian motion lies in the mean- 
ing of the coefficients, here related to the two funda- 
mental constants of the model A and a, rather than to 
the friction coefficient and the temperature of the bath. 
The quantum Brownian motion master-equation is in fact 
given by [H El El 

~A = ~\ [H , Pt )- 7 ^ [ q ,[q,p t ]} 



-l m \p,\p,Pt\ 



jl [Q,{p,Pt}] , (11) 

where 7 is the friction coefficient and f3 the inverse tem- 
perature of the background medium; the second and 
third term at r.h.s. account for diffusion, with coeffi- 
cients proportional to the squared thermal wavelength 
Ax 2 h = /3S 2 /4M and the squared thermal momentum 
spread Ap 2 h — Mj (3 satisfying the minimum uncertainty 
relation Ap th Ax th — H/2, while the last is due to fric- 
tion and ensures that the energy of the test particle 
asymptotically goes to the equipartition value depend- 
ing only on the temperature of the bath. Note that 
in the quantum description friction, which accounts for 
the finite energy growth, is of necessity related to dif- 
fusion in order to preserve the Heisenberg uncertainty 
relation [Til I20I l2l| . A fundamental result, in order 
to understand how Eq. i|l(J|) and therefore the striking 
similarity with quantum Brownian motion appears, is 
Holevo's characterization of translation-covariant gener- 
ators of quantum-dynamical semigroups |22j| . according 
to which further important restrictions can be put on 
the operators appearing in the general so-called Lindblad 
structure, once symmetry under translations is taken into 
account. In fact, according to Holevo's result, if the gen- 
erator of the dynamics C is translation-covariant, i.e., 
commutes with the action of the unitary representation 
of translations U(a) — exp[— (i/H)ap] 

C[U(a)pW(a)] = U(a)C[p]W(a) (12) 

for all real a, then its general structure, given that q 
appears at most quadratically, is the following 



-[H(p),p}+C G [p}, 



(13) 



with H(p) a self-adjoint operator only depending on the 
momentum operator p and Cg (where G stands for Gaus- 
sian) is given by 



£>g\p\ = - ^ [aoq + H ett (q,p),p] 



KpK f - - {K^K,p} 



(14) 



with 



If 



0, then one can set ipt equal to any fixed unitary vector. 



K 

H alt (q,p) 



aiq + Li{p) 
^Ui(gLi(p) 



a ,ai £ 

L\{p)q) 
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and Li(p) a generally complex function of its argument. 
The requirement of translational invariance is a natu- 
ral and compelling one for dynamical reduction models, 
since the modification of quantum mechanics by a uni- 
versal noise should by no way break the homogeneity of 
space, introducing some preferred space location. The 
restriction to mappings at most quadratic in the position 
operator q is linked to the fact that we are looking for 
a generalization of the most simple dynamical reduction 
model where A = q and the associated master-equation 
is given by 



system coupled to some reservoir, that in the commuta- 
tor at r.h.s. of l|17f) an operator appears which is the sum 
of the free Hamiltonian and some other self-adjoint oper- 
ator, this other contribution disappearing together with 
the rest of the Lindblad form when the coupling vanishes, 
as it correctly happens in Eq. (| 1 fc>|) if the fundamental 
constant A is set to zero. The general structure (|17fl can- 
not be thought of as being made up of two distinct parts, 
since the Lindblad characterization pertains to the struc- 
ture as a whole. Note however that the free Hamiltonian 
can still be put into evidence in according to 



dt 



Pt 



[Ho, ft] 



A 



[<?> [<Z> Pt}} 



(15) 



often considered in the literature (see e.g. Q| and ref- 
erences therein) even though leading to a steady energy 
increase. In view of Eq. I|14|l the most straightforward ex- 
tension of Eq. (|15|) including a friction term proportional 
to velocity is obtained setting eto = 0, i.e., no external 
constant force since we are considering the modification 
to Schrodinger dynamics for a free particle, a\ — VA and 
Li(p) = i\f\(a/h)p, thus directly obtaining the Ansatz 
given in Eq. (£1J. With this choice of functions and pa- 
rameters Eq. I|13f) gives 



1 

+A 



Aa 

H + — {q,p\,pt 



(16) 



with A = q+i (a/h) p, as in Eq. Q, which is immediately 
seen to be equivalent to Eq. I|10(l , thus giving the rationale 
for our choice for the operator A. 

Note that looking at Eq. (|l(j|) one might erroneously 
be led to think that the modification to Schrodinger dy- 
namics amounts to a change in the Hamiltonian plus a 
correction due to a universal noise given by a mapping in 
the Lindblad form with a single so-called Lindblad oper- 
ator A. This standpoint, implicitly adopted in |14| . and 
which has often led to confusion 23] , is actually mislead- 
ing. The result by Lindblad, which is strictly speaking 
only valid when the generator is bounded, but as shown 
by Lindblad himself [2(j and by the quoted results of 
Holevo also holds for the case at hand, states that the 
generator of a completely positive quantum-dynamical 
semigroup has the form 



dipt 
A 



H dt + V\(A-r t )dWt 



y ^A - 2r t A + r 2 + - (A 2 - A^ 2 )]<U 



4>t 



(18) 



even though in this equivalent expression the martingale 
structure is less evident. 



V. SINGLE-GAUSSIAN SOLUTION 

Gaussian wavefunctions are very special and they are 
often used to represent typical physical situations; we 
now show that, as for the standard Schrodinger equation, 
our stochastic equation preserves the form of Gaussian 
wavefunctions and, at the same time, we analyze their 
evolution in time. Let us then consider the following 
class of wavefunctions: 



cf> t (x) = exp \—at(x - x t ) 2 + ik t x + 7*1 



(19) 



where at and 7* are complex functions of time, while xt 
and kt are real. By following the procedure outlined in 
ref. [j, one can show that the above parameters obey 
the following stochastic differential equations 4 : 



da* 



dx t 



dkt 



2ih 



m 



— k t dt + VX 
m 



4Aa at 
1 



A 



dt. 



2a? 



dW f , 



-2\ak t dt- Vx^-dWt- 



(20) 
(21) 
(22) 



We have omitted the equation for j t since the real part 
can be absorbed into the normalization factor, while the 
imaginary part represents an irrelevant global phase. 



dt 



Pt 



(17) 



where the self-adjoint operator H is not necessarily the 
free Hamiltonian of the system, giving its dynamics when 
it is not coupled to the environment or some noise source. 
On the contrary it usually happens, e.g. when the Lind- 
blad structure appears in the reduced description of a 



A. The time evolution of at 

The parameter at is particularly important since it is 
related to the spread of the wavefunction (|19|) in position 



4 The superscripts "R" and "I" denote, respectively, the real and 
imaginary parts of the corresponding quantities. 
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and momentum, according to the following expressions: 



B. The spread in position and momentum 



a p (t) ee V(P 2 > - tp) 2 = nj {af)2 + ial)2 , (23) 



Eq. I|20(l for a t can be easily solved; one gets: 



at 



A + iB tanh ( — Bt + k 
m 



(24) 



with: 



.4 



„ Xam _ 



4A 2 a 2 rn 2 2Am 

f < ; (25) 



r, 2 



the constant /c sets the initial condition arj. 

After a tedious calculation, one can write explicitly the 
time evolution of the real and imaginary parts 5 of at- 

R mu) sin 9 sinh^it + ipi ) + cos 9 sin(c<j2i + ip 2 ) 



2V2h cosh(wii + ipi) + cos(u 2 t + ip 2 ) 



(26) 



2y/2h 



cos 9 sinh^ii + ipi) — sin# sm(uj2t + ^2) 



cosh(cjii + ipi) + cos(ti>2£ + ^2) 



2V2Xc 



(27) 



where we have introduced the following two frequencies: 

V2 W COS0 W2 = U! sin 6>, (28) 
the frequency co and the angle 9 being defined as follows: 



A^ 2 



A U 



2A ao?7io 



1CT 5 sec" 1 , (29) 
(30) 



7T 

4 ; 



note that, due to the specific dependence of both A and a 
on m as given by Eq. © , both to and 9 are independent 
of the mass of the particle, and thus are two constants 
of the model. Note also that — as it is easy to prove — 
if Oq > 0, then af > for any subsequent time t: this 
implies that Gaussian solutions do not diverge in time. 



According to (|23|l . the spread in position of the Gaus- 
sian wavefunction Ijl9|) evolves in time as follows 6 : 



/ h cosh(wii + ipi) + cos(oj2^ + ^2) 

\/2mu! sin6>sinh(cjit + ipi) + cos0sin(w2t + ^2) 

(31) 

Here we can see one of the effects of the reduction mecha- 
nism: while in the standard quantum case the spread (in 
position) of the wavefunction of a free quantum particle 
increases in time, diverging for t — > 00, the spread ac- 
cording to our model reaches the asymptotic finite value: 



<7 q (co) 



y/2>, 



muj sin f 




The asymptotic spread decreases for increasing values of 
the mass of the particle, this property entailing that, as 
we shall discuss in more detail in Sec. lVIll wavefunctions 
of macroscopic objects are almost always very well local- 
ized in space, so well that they practically behave like 
point-like particles. 

The time evolution for a p (t) can also be obtained from 
(|23|l and, as it happens for the spread in position, also 
the spread in momentum asymptotically reaches a finite 
value, which is: 



CTp(oo) 



' hmuj sin 2 9 + (cos 9 — k) 2 



2V2 



lCr 19 W — 



sin 9 
Kg m 



where: 



10 



-14 



to 



(33) 



(34) 



To conclude the section, we compute the product of 
the two asymptotic spreads: 



er ■ ov 



(cos 9 — k) 5 



1 



(35) 



which is almost equal to hj v2, with our choice JJJ and 
© for the parameters. Note however that for any choice 
of Aq and oq Heisenberg uncertainty relations are ful- 
filled. 

In accordance with |24j . any Gaussian solution hav- 
ing this asymptotic values for the spread in position and 



6 Note that the evolution of u q (t) (and also of u p (t)) is determin- 
5 The two parameters ipi and </?2 are functions of the initial con- istic and depends on the noise Wt only indirectly, through the 

ditions. constant A. 
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momentum will be called a "stationary solution" of Eq. 
JTJ . Of course, the term "stationary" does not mean that 
such wavefunctions do not evolve in time; as a matter of 
fact (see the following discussion) they always undergo a 
random motion both in position and momentum which 
never stops. The term "stationary" refers only to the 
shape of the wavefunction: stationary solutions are spe- 
cial wavefunction which are Gaussian and with a fixed 
spread in position and momentum, given by Eqs. 1321) 
and J2SJl. 



C. The mean in position and momentum 

The mean (q)t in position of the wavefunction and the 
mean (p) t in momentum satisfy the following stochastic 
differential equations which can be derived from Eqs. H21|) 
and iJUJl: 



d(q)t = — (p)tdt + ^ 
m 

d(p) t = -2Xa(p) t dt-VXh^dW, 



2a? 



dW t , 



Their average values evolve as follows: 



dt 



E[{p)t 



E[<P>t], 
-2AaE[(p) t ] 



(36) 
(37) 

(38) 
(39) 



The first equation reproduces the classical relation be- 
tween position and momentum of a particle while the 
second equation predicts that the momentum decays ex- 
ponentially in time: 

-2Xat 



E[<p)t 



Woe" 



with 



2Xa — 2AoQ!o 



1(T 20 sec" 1 , 



(40) 



(41) 



which is an extremely slow decay rate, not depending on 
the mass of the particle. 

For completeness we consider also the covariance ma- 
trix 



C{t) 



E 



(q)t-E[(q)t] 
(P>t-E[<p) t ] 



(q)t-n(q)t] 

(p>t-E[{p) t ] 



C q l(t) C qp (t) 
C pq (t) C p 2{t) 

whose coefficients satisfy the following equations: 

2 



d_ 

dt 
d 
di 



C q 2(t) 

Cqp(t) 



m 
1 



C qp (t) + X 



1 

2a? 



(42) 



C p 2(t)-2\aC qp {t)-Xh-t 



2a? 



a 



^-C p 2(t) = -4XaC p 2(t) + Xh 2 (4 



(43) 
(44) 



In Sec. IVH Bl we will discuss the physical implications of 
the above equations in connection with the dynamics of 
macroscopic objects. 



VI. ASYMPTOTIC BEHAVIOR OF THE 
GENERAL SOLUTION 

In the previous section we have seen that any Gaussian 
solution converges towards a stationary solution i.e., to- 
wards a Gaussian wavefunction with a fixed finite value 
both for the spread in position and momentum, given 
by Eqs. l-'il'l) and (|33|l . In this section we prove that the 
spread <J q {t) of any wavefunction converges with proba- 
bility one towards a q : this means that any initial wave- 
function converges to a localized solution; for the proof 
we will follow the same strategy of Ref. [lj] . 



A. The reduction process 

It is easy to see that a Gaussian stationary solution is 
an eigenstate of the operator: 



O = p — 2iha QO q, 



(45) 



where 



[sin0-i(cos0-«)]; (46) 



°° 2V2H 

the proof basically consists in showing that the variance 

a 2 (t) = (MOt-<Ot)][O-<O>]W0 ( 47 ) 

of the operator O vanishes for t — > oo. 

The first step is to re-write Oq (t) in terms of the vari- 
ances associated to the operators q and p: 

* 2 o(t) = <r 2 p (t) + ^a 2 q (t)-2^al p (t)-^, (48) 

a q ° q ^ a q 

where we have defined: 

< P (t) = \ m\[q-{q)][p-{p)m + 

(Mp-(p)][q-(q)m), (49) 

so that for a Gaussian wavefunction such as (|19|l 



2a? 



(50) 



and cr„ p corresponds to the value of cr qtP (t) when the 
state \ipt) is a stationary Gaussian solution. 
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After a rather long calculation (see Appendix A for the 
details), it is possible to show that: 



dt 



E[4(i)] = 



-4AE 



_2 2^,-4 (rfW <p(t) 



1,P^ 
—2 
U 1,V 



< 0. 



(51) 



Since <Jq (t) is by definition a positive quantity, the above 
equation is consistent if and only if the r.h.s vanishes for 
any lo S f2, with the possible exception of a subset of 
measure 0. This is particular implies both that a 2 (t) — > 
a.s. and that cr q (t) — ► ~d q a.s., which is the desired result, 
i.e., the wavefunction converges to a localized solution. 



B. The localization probability 

Once proved that Eq. Q with the choice (0J for the 
operator A induces the localization of the wavefunction 



in space, it becomes natural to analyze the probability 
for a localization to occur within a given interval of the 
real axis. Such an analysis can be developed along the 
same line of Ref . £| . 

Let us consider the probability measure: 



Mt(A) = 



(52) 



defined on the Borel sigma-algebra B(M.) of the real axis, 
where Pa is the projection operator associated to the 
Borel subset A of M; such a measure is identified by the 
density p t (x) = E P [\ip t (x)\ 2 }: 



Mt(A) 



p t (x) dx. 



(53) 



A 



The density pt(x) corresponds to the diagonal element 
(x\pt\x) of the statistical operator p t = Ep [|i/) t )('0 t |], so- 
lution of the master-equation I|1U|) . In Appendix [B] we 
show how the general solution of this master-equation 
in the position representation can be obtained; the fi- 
nal expression, as a function of the solution of the free 
Schrodinger equation pf (i.e., with A = a = 0), is: 



\Q1\Pt\q2) 



dk / dye~^^ k yF[k, qi ~q 2l t] 



x ( y + 



1 



2m \ 2Xat J 



PI 



gi + 92 _ qi - g 2 c - 2 \ at 



kt 
2m 



2\at J I 



(54) 



with 

F[k, x, t] = exp ■ 

and 



Aa 2 

2h 2 



kH 



8aT 2 



-2\at 



r t 



2m\c 



k Ki(t) +2x [x 



-2\at 



2m\a 



k)K 2 (t)+x 2 K 3 {t) 



Ki{t) = T 2 + 2T t - AXat, 

K 2 {t) = e- iXat + 4Xate~ 2Xat - 1, 

K 3 (t) = -4\ate- 4Xat —T 2 + 2T t e- 2Xa \ 

where we have defined T t — 1 — e~ 2Xat . Taking the diagonal matrix elements one has: 



(55) 



(56) 



1 



Pi I 



2nh 



( x ) = 7TT dk d V ^ {l,h)kv 0, t] ( y + x + — 1 



kt 



2m 



and 



F[k,0,t] =exp<j~^fc 2 t 



1 



2Xat J 
■k 2 Ki(t) 



y + x 



kt 
2m 



2Xat J 



(57) 



(58) 



According to the values (0) and JSJ) for Aq and a>o, the above expressions can be expanded for small 7 Xat, leading 



7 This means that we are considering only times t <C (Aa) 1 ~ 
10 20 s. 
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to: 



Pt(x) 



1 

2nh 



dk 



dye-W) 1 * exp 



Xk 2 
6m- 



\a 2 k 2 
2h 2 



y + x + 



Xakt 2 
2m 



y + x- 



Xakt 2 
2m 



(59) 



We now focus on the case of a macroscopic object (let us say m > 1 g); one can further approximate the above 
expression by noting that for such values of m the exponential factor appearing in Eq. H59[l damps all matrix elements 
such that the term Xakt 2 /2m is not vanishingly small; e.g. when Xakt 2 /2m > 10~ 15 m, then the second exponential 
in the above equation is much smaller than e~ 10 (™ 1 / k k) We can then neglect the two terms in the matrix elements 
and perform the integration over k, and we get: 



Pt(x) 



dye~ 0ty pf{x + y) 



(60) 



with 



2h 2 



i 3 



10 



43 



> 10' 



(™) for:t>10- n 



kg J V t J 



s. 



„ , m \ /sec 

65 / 1 ( 



kgj V t > 



\ for: t < 1Q- 11 s 



(61) 



The exponent in (|60|l is extremely peaked with respect 
to the typical values the probability density pf(x) asso- 
ciated to the wavefunction of a macroscopic object takes, 
so that with very high accuracy we have: 

p t (x) ~ pf(x). (62) 

As discussed in ref. 0, the probability measure /xt(A) 
which we have shown to be extremely close to the quan- 
tum probability obtained from the free Schrodinger equa- 
tion can be interpreted as a probability measure for the 
collapse of the wavefunction of the macroscopic object 
within A, when A corresponds to an interval of the real 
axis of width greater or equal e.g. to 10 -5 cm |4j. 

To conclude, the previous analysis shows that under 
the above listed conditions the probability for the wave- 
function of a macro-object to be localized within an inter- 
val of the real axis is almost equal to the corresponding 
quantum probability as given by the Born rule. 



VII. DYNAMICS OF MICROSCOPIC AND 
MACROSCOPIC SYSTEMS 

In this section we discuss how our reduction model is 
related both to quantum and to classical mechanics. Our 
aim is to show that at the microscopic level the physi- 
cal predictions of the model are almost identical to stan- 
dard quantum predictions and that, at the same time, the 
model with high accuracy reproduces classical mechanics 
at the macro-level. 



A. Micro—systems: comparison with standard 

quantum mechanics 

Microscopic system cannot be directly observed, and 
their properties can be analyzed only by resorting to suit- 
able measurement procedures. All physical predictions 
of our model, concerning the outcome of measurements, 
have the form Ep[(-0 t |S'|-0 t )] where S is a suitable self- 
adjoint operator, typically a projection operator and it 
is easy to show that Ep[(ip t \S\tpt}} = Tr[Sp t ], where the 
statistical operator p t = Ep[|^)('i/; t |] satisfies Eq. l(TU|l . 
Accordingly, as already discussed in Sec. lIVI the testable 
effects of the stochastic process on the wavefunction are 
similar to the effect induced by a quantum environment 
on the particle, when both friction and diffusion are taken 
into account. 

With our choice J7J and (JSJ for Ao and ao, the testable 
effects of the stochastic process are of the same order 
of magnitude of those induced by the interaction of the 
system with particles and radiation of the intergalactic 
space 25]: such effects are very small and masked by 
most other sources of decoherence, so that they can be 
tested only by resorting to sophisticated experiments j2|| 
|27| . This implies that the physical predictions of our 
model are very close to standard quantum mechanical 
predictions. 

B. Macro— objects: comparison with classical 

mechanics 

A macroscopic object is made of elementary con- 
stituents strongly interacting among each other and, ac- 
cording to our model, its dynamics is governed by the 
following stochastic differential equation, which is the 
straightforward generalization of Eq. to a system 
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of N particles: 



i 



N 



H TOT dt + J2 V\i {A n - r nt ) dW t n 



N 



<MW); 
(63) 



the symbol {x} = xx, %2, ■ ■ ■ > %n represents the N spatial 
coordinates of the configuration space of the composite 
system; W} , , . . . Wf* are N independent Wiener pro- 
cesses; r nt — (tjj t \[A n +Al l )\ipt}/2 and the localization op- 
erators A n are given by expressions @, with q replaced 
by q n , the position operator of the n— th particle, andp re- 
placed by p n , the corresponding operator. Furthermore: 



and A CM and a CM defined by Eqs. I©, with to equal to 
the total mass M of the composite system. Note that 
the separation of the center-of-mass motion from the rel- 
ative motion, for the non-Schrodinger terms of Eq. i|63[l . 
is possible because of the specific dependence of the pa- 
rameters A n and a n on the masses m n of the particles as 
given by Eq. 

According to Eq. (|67|l the center of mass behaves like 
a particle whose dynamics, in the free case, has been 
discussed in detail in Sees. IVland lVll we now show how 
the large numerical value for M, typical of macroscopic 
objects, affects the time evolution of the center-of-mass 
wavefunction. 

1 Collapse rate. According to Eq. I]5ip. and as- 
suming that the wavefunction is not already localized in 
space, i.e., that <r q (t) 3> ~5 q , one has: 



and 



to 



to 



(64) 



dt 



n*o(t)] 



> 



0~Z 



2A cj 2 sin 2 9 
to 



m 3 a 4 Jt), (71) 



N , 

2 

n=l 



Hq being the standard quantum Hamiltonian for the 
composite system. 

As custom, we separate the motion of the center of 
mass from the internal motion. To this end let: 



1 N ( 

71=1 \ 



{q n ,Pn}, (65) with 2X uj 2 sin 2 6/m ^ 10 15 Kg" 1 m" 2 sec" 3 . In the 
microscopic case the r.h.s. of Eq. I|71|) is in general small 
and negligible as we expect it to be, since the reduction 
mechanism must be ineffective on microscopic systems; 
in the macroscopic case instead it rapidly becomes big, 
due to the large value of the mass m, ensuring that 
any initial wavefunction rapidly converges towards a 
localized solution. We can then assume that, possibly 
after a very short transient period, the wavefunction 
describing the motion of the center of mass of any 
macro-object is practically localized in space. 



N 

M — m Tl 

n=l 



(66) 



be the position operator associated to the center-of-mass 
coordinate R, and q n the position operators associated 
to the internal coordinates x n = x n — R (n = 1, . . . , N); 
let also P and p n be the corresponding momentum oper- 
ators. Then, if H® OT can be written as the sum of a term 
if° M associated to the center of mass and a term i?° EL 
associated to the internal motion, it is easy to prove that 
the dynamics of the two types of degrees of freedom de- 
couple; in particular the equation for the center of mass 
— the only one we consider here — becomes: 



2 Behavior of the stationary solution for macro- 
scopic objects. We now discuss the time evolution 
of a typical localized wavefunction, i.e., a Gaussian 
stationary solution. Eqs. (|36fl and (|37|l imply that the 
two maxima in position and momentum of such wave- 
functions fluctuate around their mean values; we now 
show that in the macroscopic regime these fluctuations 
are extremely small. 

As a matter of fact Eqs. (|42|l to H44|) imply for a sta- 
tionary solution: 



dMR) = 

A C M 



2 ("CM^CM 



2r CM|t A CM + r~ 



dt 



with: 



A, 



c. = 2<M4 



+ 1- 



CM 

a C M 



{Q,P}- 



+ ^cm]|V>*>, 



P. 



t ,t)dW t 
MR), (67) 

(68) 
(69) 
(70) 



C,2(t) 



xft 



hf 



2Xa 2 m V sin 



cos v — K 



(1 



-2Xat\ 



C P 2 ( t ) = TZ 



16A 2 a 3 TO 2 

h 2 



COS u — K 



COS V — K 

4a \ sin 9 
where I is defined as follows 

h 



sm( 

2 



(1- 



1-e 



-4\at 



e = - 



2\am V sin 



cos u — K 



V2h 



tow sin ( 



) (72) 



(73) 



(74) 
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Since the exponential factors decay very slowly (we re- 
member that Xa ~ 10~ 20 sec -1 ) it is physically signifi- 
cant to consider only the linear term of their Taylor ex- 
pansion; one then gets: 



C q *{t) 

C p 2 (t) 



Ah 2 \ t 

moLu 2 m 

h 2 x 



10 



-33 



Kg 



t 



m 



m 



-mt 



10 



-43 



m \ 



sec^ 



(75) 



Kg 2 nr 
Kg J \sec ) sec 2 



(76) 



which are very small quantities, when m is the mass of a 
macro-object. Accordingly, in the macroscopic case the 
actual value of the two peaks in position and momen- 
tum of a stationary Gaussian solution are very close to 
their average values which, as we have seen, evolve in 
time according to Newton's laws for a free particle mov- 
ing in a (very weakly) dissipative medium. This proves 
that with high accuracy a stationary solutions for the 
center of mass of a macro-system practically behave like 
a point moving in space according to the laws of classical 
mechanics. 



VIII. TIME EVOLUTION OF THE MEAN 
ENERGY 

We now discuss one of the main purposes of our work, 
i.e., we show that within our collapse model the energy 
of an isolated system does not increase with a constant 
(even if negligible) rate, but reaches an asymptotic finite 
value. Indeed, this result is entailed by Eq. (|10|l for the 
statistical operator, which motivated our choice for the 
localization operator A. 

As a matter of fact, the mean value K[(H ) t ] of the 
energy satisfies the following equation: 



d \h 2 

j t n(H )t] = ^ - 4XaE[(H ) t ] 



whose solution is: 

E[(ff >t] = I Eu 



8ma 



-AXat 



8ma 



(77) 



(78) 



As we see, the mean energy of an isolated system changes 
in time and thus is not conserved; anyway it does not 
diverge for t — ► +oo, but reaches the asymptotic finite 
value: 



of change being equal to AXa — 4Ao<2o — 10~ 20 sec -1 
which means that, with very high accuracy, the mean 
energy remains constant for very long times. 

The above equations imply that the stochastic pro- 
cess acts like a dissipative medium which, due to fric- 
tion, slowly thermalizes all systems to the temperature 
T by absorbing or transferring energy to them depend- 
ing on their initial state. Note that according to Eq. (HQ)) 
a very "cold" medium is enough to guarantee the local- 
ization in space of the wavefunctions of macroscopic ob- 
jects. Note also that one recovers a GRW-type equation 
by setting a — > 0, which corresponds to the high tem- 
perature limit T —t +oo. This implies that the reason 
why in the original GRW reduction model the en- 
ergy increases and eventually diverges is that the noise 
acts like an infinite-temperature stochastic medium. Our 
model also shows that the increase of the mean energy of 
a quantum system subject to spontaneous localizations is 
not an intrinsic feature of space-collapse models (indeed, 
according to our model, the mean energy decreases for 
typical quantum systems) and it can be (partly) avoided 
with a suitable choice of the localizations operators. 

The above discussion suggests that the model can be 
further developed by promoting Wt to a real physical 
medium with its own equations of motions, having a 
stochastic behavior to which a temperature T can be 
associated and such that, with good accuracy, can be 
treated like a Wiener process. This would imply not only 
that the medium acts on the wavefunction by changing its 
state according to Eq. QJ, but also that the wavefunction 
acts back on the medium according to equations which 
still have to be studied. The above suggestion opens the 
way to the possibility that by taking into account the 
energy of both the quantum system and the stochastic 
medium one can restore perfect energy conservation not 
only on the average but also for single realizations of the 
stochastic process. A similar proposal has been consid- 
ered by P. Pearle [28| and by S. Adler 0. 

These ideas will be subject of future research; we con- 
clude by noting that, whatever its nature can be, the 
stochastic medium cannot be quantum in the usual sense 
since its coupling to the particle is not a standard cou- 
pling between two quantum systems: Eq. Q, in fact, is 
not the standard Schrodinger equation with a stochastic 
potential. 



IX. CONCLUSIONS 



lim E[(H )t] 

t — > + OG 



8ma 

corresponding to a temperature 

n 2 



T 



8mo«o ' 



K, 



(79) 



(80) 



which is independent of the mass of the particle. Note 
that the time evolution of E[(Ho) t ] is very slow, the rate 



We have presented and analyzed a collapse-model 
which preserves the standard quantum mechanical pre- 
dictions and reproduces classical mechanics at the macro- 
scopic level, at the same time avoiding the infinite growth 
of energy of the system, a criticized feature of the 
space collapse models appeared in the literature. This 
has been obtained by drawing on an analogy with quan- 
tum Brownian motion, where friction effects, directly re- 
lated to diffusion for the preservation of Heisenberg un- 
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certainty relation, guarantee that the energy of the test 
particle reaches a finite value depending on the parame- 
ters of the model. The model is also characterized by the 
fact that the related master-equation complies with the 
general structure of translation-covariant generators of 
quantum-dynamical semigroups obtained by Holevo |22j | , 
so that symmetry under translation is correctly kept into 
account, as compulsory for dynamical reduction models, 
which should not introduce any preferred space location. 

Needless to say the exploited analogy with the quan- 
tum Brownian motion master-equation, typically used for 
the description of dissipation and decoherence, should 
by no way induce confusion on the different nature of 
the two issues of decoherence and of the measurement 
or macro-objectification problem in quantum mechan- 
ics: as it has been stressed also in recent publications 
|25l l29| decoherence does not provide a solution to the 
measurement problem. This important conceptual differ- 
ence notwithstanding, dynamical reduction models and 
models of environmental decoherence both impinge on 
the same mathematical inventory, typically used in the 
theory of open quantum system 30], so that results ob- 
tained in the one field can often be fruitfully exploited in 
the other. 

An extension of this approach might be pursued in or- 
der to cope with the infinite energy growth also in the 
original GRW model of dynamical reductions Q, build- 
ing on the analogy with the quantum linear Boltzmann 
equation plj . 
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APPENDIX A: THE MEAN RATE OF CHANGE 
OF THE VARIANCE a 

In order to derive Eq. I|51|l , we first compute the aver- 
age value of the time derivative of the second moments of 
the operators q and p and of their symmetrized correla- 
tion. Using the evolution equation , we obtain through 
Ito calculus: 



2E\a 2 (t)} 

- Kit). = [ ™ K n -4AE[^(t)] + 4aAE[^ft)], 

(Al) 



d_ 

It 



- E[a 2 (t)} = -4AE[<4,(t)] - 4aAE[o*(t)] + \h 2 , 

(A2) 

d o E[cr^(i)l 

-E[< p (i)] = J-^21-AXE[al p (t)al(t)}. (A3) 



Since at is constant for a stationary solution, the station- 
ary values of variances and correlation are such that the 
right-hand sides of I|A1|I to (|A3|) vanish: 



2 

^■-2Xat + 2aXa'i = 0, 
m q q 



\ —4 i 2c ^ ~ 2 

_ _ h 2 

<Tq,p + a(J l ~ -J = °> 

^ - AXa 2 , Sal = 0. 



Moreover Eq. (35J) can be rewritten as: 



2 2 4 



(A4) 
(A5) 
(A6) 

(A7) 



For convenience we define the quantity X, Yand Z by 
the equations: 

a 2 (t) = a 2 q (l + X), 
a 2 p (t) = o*(l + Y), 
°l P (t) = 5t P (l + Z), 
so that Eq. (|48|l becomes 

2at 



a 2 Q {t) = a 2 p {X + Y)-^fZ, 



(A8) 



where we also made use of the Eq. <|A7ll . 

Using Eqs. (53J to JSHJ we obtain for the average 
value of the time derivative of <Jq : 



dt 



E ko(*)] = WiE[X] + w 2 E[Y] + w 3 E[Z] 



- 4A (a 2 q a 2 p E[X 2 } + a% p E[Z 2 ] - 2a\^E[XZ\) (A9) 



with 



-8A ( alal aaZ - a 



^2 —4 

YT — CI 



w 2 = -2 ( 2a\a 2 ' °* 



P ' — 2 I ' 



w 3 



2 2 

' m a 2 q 



From Eqs. HA7jl and l|A5|l we get for w% 



Wl 



-4Xa 2 q a 2 p , 



while using Eq. (|A4() one can prove that 



W 2 = Wl. 



(A10) 
(All) 
(A12) 

(A13) 
(A14) 



Finally, from Eq. I|A6(I . we find that w 3 and w\ are re- 
lated as: 



2a 



w 3 



r=4 
^f2— 2 



(A15) 
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so that we have 



Wl X+w 2 Y+w 3 Z = Wl [x + Y-^Z) = -4A<x^($, ecomes 



corresponding to the characteristic function associated 
to the Wigner function. The master-equation IjBlfl thus 



(A16) 

where the last equality is obtained through Eqs. (|A8(I 
and (|A13|I . Finally, the above result together with Eq. 
(I A 71) allow us to write 



j t nrt(t)] = -4AE 



which is just Eq. 



ay o {t) + al p {X-Zf+'^X 2 



n 2 

4 

(A17) 



° ~ n \ k ® ~ n \ ^ 2~ n \ 
—p t (k,x) = ——pt(k,x) - -x pt(k,x) 

at m ox I 

\o? d 
- -^k 2 p t (k,x) - 2Xax— p t (k,x), (B6) 



while HB2I) and HB4I) take the form 



APPENDIX B: GENERAL SOLUTION OF THE 
MASTER-EQUATION 

We now show how to obtain the general solution of the 
master-equation 

j t Pt = ~\ [H ,pt]-- [q,[q, Pt ]] 



pt(k,x) 
and 
p t (k,x) 



dk dx G{k, x, t; k , x , 0)p o {k o , x Q ) (B7) 



dkodxQ j drds G(k, x, t; ko, xq, 0) 

x G o {k o ,x o ,0;r,s,t)pf{r,s) (B8) 



Xa 2 
2ft 5 " 



\p,\p,p t }] - i^- [q,{p,Pt}}, (Bl) 



respectively, where pf(r,s) is again the solution of the 
free Schrodingcr equation, Go is simply given by 

G (k, x, t; ko, x ,0) = S(k— k )S ( x — x + —t) , (B9) 



G{k, x, t; ko, xq, 0) — ^> S(k — ko)S(x — xq). (BIO) 



given in 1)10(1 . which has the same form of the quantum 

Brownian motion master-equation, partially following an and G satisfies the following initial condition: 
appendix of |25j . In particular we want to express the 
general solution in the position representation {q\\pt\q2) 
as a function of the solution of the pure Schrodinger 
equation (qi\pf 1 52} , in order to point out the corrections 
to the position probability density due to the non-linear 
stochastic modification of the Schrodinger dynamics. 

As a first step we want to express the solution of I|B1|) 
as a function of the generic initial condition po according 
to 



Ea. (|B8l) can be brought back to (|B4() by exploiting the 
inverse relation of Eq. I|B5|) . i.e., 



{Ql\PtW) 



1 

2^h 



dke'^ k ( S1 ^ L )p t (k,q 1 ~q 2 ). (Bll) 



(Ql\Pt\Q2) = J dq w dq2oG(qi,q2,t;q 10 ,q2o,0)(q w \po\q2o}, 

(B2) 

where G is the Green function solution of the partial 
differential equation associated to i|Bl|) in the position 
representation satisfying the following initial condition 

G(qi,q 2 , t; q w , q 20 , 0) — S(qi - qw)S(q 2 - q 20 )- (B3) 

Once G is known one may immediately express the 
general solution as a function of the solution of the 
Schrodinger equation by means of the free propagator 
Go 



\q1\Pt\q2) 



The key observation in order to determine G is the fact 
that Eq.jBlJ) preserves Gaussian states, so that given the 
Gaussian Ansatz: 

pt{k, x) = exp {— c\k 2 — c 2 kx — c$x 2 — ic^k — ic$x — c§} 

(B12) 

one easily obtains the evolved state by expressing the 
coefficients at time t, Ci(t), as a function of the initial 
state characterized by the values of the coefficient at time 
zero. Considering an initial Gaussian state pQ° x "' ev (k,x) 
which in the limit of small e and 77 approximates the Dirac 
delta, according to 



pl m {k,x) 



dqwdq 20 j dudvG(qi,q 2 ,t;q w ,q 20 ,0) 
x G a (q w ,q 20 ,0;u,v,t)(u\pf\v). (B4) 



For calculational purposes it is however convenient to 
consider the quantity 



£,7)->0 



► 5(k — ko)S(x — xq) 



one has 



pt(k,x) = Tr (pte* 



(kq+xp) 



(B5) 



(B13) 



(B14) 



Coming back to the Ansatz i|B12J) , the coefficients satisfy 
the equations 



G(k, x, t; ko, Xq, 0) = lim p koX "' er, (k,x) 

e,ri—>0 
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ci(t) 
<*(*) 

Ci(t) 



c 2 (t) | Aa 2 
m 2H 2 



2c 3 (t) 



2Xac 2 (t) 



X 
2 

<%(*) 



4Xac 3 (t) 



(B15) 



with solutions 



ci(t) = ci(0)+ca(0) 



2mAa 



cs(*) = -2Aac 5 (t) 
-c 3 (Q) 



, , - QO * , (r 2 + 2r t - AXat) + 

im 2 X 2 a 2 32m 2 X 2 a 3 2h 2 



c 2 W = c 2 (0) e - 2AQt + c 3 (0) 



p -2Aai 



+ 



Xa 8mXa 2 



1 



c 3 (t) = -+^(0)- 
c 4 (t) = c 4 (0) + c 5 (0) 



1 



-4Acrf 



(B16) 



2mAa 



c 5 (*) = c 5 (0)e 



-2Aat 



,-2Aat 



with r* = 1 



one has, exploiting (|B16I 



and cq simply a constant. For the choice 



Po ox °' ev (k,x) 



i 



i 



i 



1 



exp < — (fc — ko) > exp < (x — xq) 



(B17) 



exp <| ( /.' — An I" J> exp < 



_£tk_ 

2mXa 



- (x - xe- 2Xat ) 



Xa 1 



x exp < — -—r k t> exp < — 



2h 2 



1 



8n 



^(l- e - 4AQt ) 



4to 2 A 2 q; 2 



mXa 



(B18) 



so that taking the limit one has 

G(k, x, t; kg, xq,0) = 5(k—ko)S 

with 



Ttk 

2mXa 



{x - xe~ ZXat ) exp 



Aa 2 



k t > exp 



J -^—[xlK^t) + 2xx K 2 (t) + x 2 K 3 (t)} 



[8aT 



Kt(t) = T 2 t +2T f -AXat 
K 2 {t) = e- iXat + 4Xate- 2Xat 



- 1 



(B19) 



(B20) 



K 3 (t) = -4Xate~ 4Xat -T 2 t + 2Y t e 



-2Xat 



and therefore (|B8f) now explicitly becomes 



Pt(k,x) = exp 



Aa 2 
2h? 



k t > exp 



8aT 2 



xe 



-2\at 



T t k 
2mXa 



Kt(t)+2x ixe- 2Xat + 



J A.,; , r tk \ 

2mXa J 



K 2 (t)+x 2 K 3 (t) 



x pf ( k,xe' lXat + 



-2Xat , 1 < 



T t k 



2mXa 



2Xat\ 



Exploiting the inversion formula 1|B11|I together with the expression 

Pt(k,x) = J dye^ ky {y + ^\p t \y-^) 



(B21) 



(B22) 
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equivalent to (|B5|) one finally obtains the desired explicit expression for (|B4(1 : 

Xa 2 



2h 2 



k z t 



xexp 



{sar? 



(qi - q2)e 



~2\at 



2m\a 



k K 1 (t) + 2(q 1 -q 2 ) (?i-fla)e 



-2Xat 



2m\a 



k K 2 {t) + ( Ql - q 2 ) 2 K 3 (t) 



1 - 



2Xat J 



gl +g2 _ gl ~g2 2Aot _ 

,y 2 2 2m 



2Acrf 



)>• (B23) 



r 
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